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Abstract: The inflaton might be coupled to a gauge field through a term 
f'^{(t))Ffj_„F^^ . If / oc where a{t) is the scale factor, the perturbation 5W of the 
gauge field generates a potentially observable statistically anisotropic contribution 
to the primordial curvature perturbation during slow-roll inflation. The spectrum 
and bispectrum of this contribution have been calculated using the in-in formalism 
of quantum field theory. We give a simpler and more complete calculation using 
only the classical perturbations. The results suggest that either the entire curvature 
perturbation C, (both the statistically isotropic and anisotropic parts) is generated 
during slow-roll infiation, or else it is generated afterwards. 
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1. Introduction 

The observed primordial curvature perturbation ({x.) presumably originates in the 
early universe, and a central task of theoretical cosmology is to determine its origin 
and the evolution ({-x,t). It is usually supposed to be generated from the pertur- 
bation of one or more scalar fields, making it statistically isotropic which is con- 
sistent with observation. But it may also receive a contribution from one or more 
vector fields, making it statistically anisotropic. Observation does allow significant 
anisotropy and in particular a spectrum of the form^ 

Vc{k)=Vc{k) [l + g^iW-kf] (1.1) 

is allowed with l^f*! < 10~^. The upcoming PLANCK results will give < 10~^ 
barring a detection ||1|. 



We use hats to indicate unit vectors 
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Most schemes invoke a U{1) gauge field B^. If tlie field has the canonical kinetic 
term and no coupling to gravity, its spectrum increases with wavenumber like k"^, 
and its contribution to the spectrum of ( has the same behaviour making it almost 
certainly negligible on cosmological scales 0. One way of avoiding this may be to 
keep the canonical kinetic term but invoke a coupling to gravity given by —RB^B^'^/ 6. 
This apparently |]^, ^ generates a flat spectrum but the model has instability and 
(possibly) non-linearity [Q and it is not known [^] whether or not these spoil the 
prediction. 

In this paper we consider a different scheme, which invokes a gauge kinetic func- 
tion that depends on the slowly rolling infiaton field 0. To allow a significant 
contribution on cosmological scales one needs (at least approximately) / oc a^^ where 
a{t) is the scale factor, which may be reasonable in the context of string theory [Q. A 
well-defined contribution to C is generated during slow-roll inflation 0, || and addi- 
tional contributions may be generated afterwards [p, In this paper we give for 
the first time a master equation that includes all contributions, but our main focus is 
on the contribution generated during slow-roll inflation. Under the assumption that 
the homogeneous part of the gauge field dominates the perturbation, it has been cal- 
culated [|T^, § using the in-in formalism of quantum field theory. Using instead the 
classical formula for C, in terms of the energy density perturbation, we reproduce that 
calculation and extend it to the case where the homogeneous part is sub-dominant 
or even negligible. We also identify the assumptions made by the calculation, and 
show how to perform a more complete calculation with weaker assumptions. We end 
by discussing the significance of our result. 

While this paper was in its final stages of preparation ref. [11 appeared. The au- 
thors of this work uses the so called 5N formalism to calculate and non-gaussianity 
in the same model. Our results agree with their findings where there is an overlap. 



2. Action and field equations 

During inflation we are taking the action to be 



S 



]^MlR - ^9^09^0 - F(0) - \f{ct>)F^,F^'' 



(2.1) 



where F^j^ = 9^-6^ — dyB^j_ with i?^ a gauge field, and / is some function of the 
infiaton field 0. 

Extremizing the action in Eq. (|2.1| ) with respect to fields 0, B^ and their deriva- 
tives we obtain field equations 



[d, + d. In ^] 9> + V' + -ff'F.^F^'^ = 0; 

[d, + d^\n^]fF'^'' = 0, 



(2.2) 
(2.3) 
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where g = det{g^u), V = dV/d(j) and /' = df/d(j). We make the gauge choice 
Bq = 0, which fixes the spatial components Bi{-x,t) up to a constant. 

We are going to assume that the infiationary expansion is nearly isotropic, check- 
ing later that this is justified. Writing 



we then have 



ds'^ = -dt^ + a^{t)Sijdx'dx^, (2.4) 



+ 3H^ - a-^V^<f) + V' = -^ff'F^.F^", (2.5) 
Bi+\H + 2YiB,- a-^V^Bi = a"^2^d,Bi, (2.6) 



where H = a/a and = bi^d'^ jdx^dxK 



3. Pure slow-roll inflation and the curvature perturbation 

If the effect of the gauge field is completely negligible we have pure slow-roll inflation, 
described for instance in |jT2|, |13]. The last term of Eq. (|2.2| ) is absent, and for the 
unperturbed inflaton (/)(t) we have 

+ 3if0 + 1/'(0) = 0. (3.1) 

The potential is supposed to satisfy the flatness conditions e ^ 1 and |?7| <^ 1 where 
e = M|(y/V)^/2 and ?7 = M^V" jV . Then, more or less independently of the initial 
condition 

3i/0 ~ -V. (3.2) 

With the flatness conditions this is called the slow-roll approximation, which we use 
throughout. It gives \B.\jH ^ H and |e|//7 ^ e. Except where stated we take if 
and e to be constant. 

We write </)(x, t) = <^(t) + 5</)(x, t) and deflne 

6Mt) = J d^xe-'^-''6(p{^,t), (3.3) 

and similarly for other quantities. 

Let N{k) be the number of e-folds of slow-roll inflation after the epoch of horizon 
exit k = aH. The 'cosmological scales' on which (^k is observed have A^(/co) ~" 15 ^ 
N{k) < N{kQ), where ko = aoHo and the subscript on the right hand side denotes 
the present epoch. Also, 

15 < N{ko) < 70, (3.4) 
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where the upper bound assumes P < p after inflation, were P is the pressure and p 
the energy density. The lower bound is needed so that cosmological scales leave the 
horizon. Typical cosmologies give a value in the upper third of the range. 

Perturbing Eq. ( p.2|) to first order, and ignoring the metric perturbation (back- 
reaction) we have 

(50k + 3i/(50k + (fc/a)^50k = -l^"50k. (3.5) 
Evaluating the metric perturbation to first order, on the flat slicing of spacetime 



(such that the 3-curvature scalar vanishes) one finds [|T4|, |T3] that it is significant 
only after horizon exit, when Eq. (|3.5D becomes 



6(f) + 3/f(50k ^ -V"6(j)^ - 6H'e6(f)^. (3.6) 

The last term represents effect of the metric perturbation (back-reaction). From 
Eq. (|2.2| ) this term can be written 

a~^do6{\n v^)a^0 ^ 3H6{\n 7^)0. (3.7) 

Comparing Eqs. (|3.6|) and (p.7|) we find^ 



5(lny^) ~ ^^50k. (3. 



Regarding (50k(t) as an operator its mode function (f){k,t) satisfies Eq. 
Well before horizon exit at A; = aH Eqs. ( |2.1| ) and ( p.2|) describe a free field in flat 
spacetime. We choose 

a{t)(l){k, t) = e-'''/^^/V2k, (3.9) 
and the vacuum state corresponding to the absence of </> particles. Then we have 

(50k50k') = {27:f6'\k + k'){27ry k')V4k,t), (3.10) 

{271^ e)V^{k, t) = (l)\k,t). (3.11) 

Also, 

i^M^JM = 0, (3.12) 

k2^0k3(^0k4) = (50ki50k2)('^0k3^0k4) + 

+ (50ki50k3)('^0k4'^0k2) + ('^0ki50k4)('^0k3'^0k2), (3.13) 

and similarly for higher products. 

Well after horizon exit [k -C aH <^ 1) the phase of (f){k, t) becomes constant 
which means that 50 can be regarded as a classical perturbation with the above 
correlators (gaussian perturbation). In the initial regime 

\r]\,e<^k/aH <^1, (3.14) 



^We assume t/) < as is the case for hybrid inflation. 
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the left hand side of Eq. ( |3.5|) is still close to zero, and the mode function is 

t) ^ Mk, t) ^ (1 - zaH/k) ^ (3.15) 
This gives a slowly varying perturbation, H~^\6<pk\ <^ |^0k| with 

P,(A;,t)^pJ^ (^^y. (3.16) 



The subsequent evolution is given by Eq. ( ^.6| ) . It has two independent solutions but 
the initial slow variation will pick out the solution 

3H6<j)k - - [V" + QH^e] (50k. (3.17) 

Going to second order in the perturbations of the field and metric one finds addi- 
tional correlators of 5(j) (non-gaussianity) [0, starting with the three-point correlator 
which defines the bispectrum B^: 

{5(t>^,M^^M^^^) = (27r)353(ki + ks + k3)5^(ki, ks, kg). (3.18) 

To define the primordial curvature perturbation (^(x, t), one smoothes the metric 
and the energy-momentum tensor on a super-horizon scale and chooses the comoving 
threading and the uniform energy-density slicing of spacetime. Then 

C(x,t) = 51n(a(x,t) = 5[lna(x,t)/a(ti)] = 5[A^(x, t, ti)], (3.19) 

where a(x, t) is the local scale factor such that a comoving volume element is pro- 
portional to a^, and a{ti) is its unperturbed value which can be evaluated at any 
epoch. As we discuss in Section ^]5| it is usually enough to work to first order in 
then one can choose ti = t to get 

C(x, t) = HSt{x, t) = -H5p{x, t)/p{t), (3.20) 

where 6t is the displacement of the uniform-p slice from the flat slice and 6p is 
evaluated on the fiat slice and 

p = -3H{p + P). (3.21) 

Since p is smoothed on a super-horizon scale, Eq. ( |3.21|) applies at each location, 
which implies that C(x, t) is constant during any era when P(x, t) is a unique function 
of p(x, t). Galaxy surveys and observation of the CMB anisotropy detect the value 
of ( at the epoch with temperature T ~ 10~^ MeV; the universe is then radiation- 
dominated (P = p/3) giving ( a constant value that we denote simply by ({'x). Its 
spectrum is nearly scale- independent with [|T7| ^ 

V^{k) ~ (5 X 10-^)2 (3.22) 
n{k) - 1 = d\nV^{k)/d\nk ~ -0.029 ± 0.010. (3.23) 
The spectrum and bispectrum of any perturbation are defined as in Eqs. and ( |3l8| ). 
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The result for n{k) assumes that it has neghgible scale dependence. It also assumes 
a tensor fraction r <^ 10~^, which will soon be tested by PLANCK [p!8| .) 
For the reduced bispectrum defined by 

5 B 

where P({k) = {2tt^ /k^)V(, current observation J7[ give |/nl| ^ 100 and barring a 
detection PLANCK will give |/nl| ^ 10. For /nl to ever be observable we need 
I/nlI > 1. 

During pure slow-roll inflation, ( has a time-independent value that we denote 
by It is well- approximated by Eq. (|3.20| ) with 

p^^V(0) + i02^\/(0), (3.25) 

leading to the gaussian perturbation ( = = —Hdcp/cj) and 

1 fH 



where e and H are evaluated at horizon exit. 



Going to higher order in the perturbation of the metric and of 50, one finds [115 



a non-zero B^. Then, using the 5N formula to go beyond the first-order formula 
(PD, one can use Eqs. [ WM . O, and ( gl3D to find |/nl| ~ {V(^jV(^flQ~^ 
which is too small to observe [IlSl |l9l.^ 



4. Background (unperturbed) universe 



Recasting Eqs. ( p.SP and ( p.6P in terms of W with components Wi = fBi/a and 
dropping gradient terms, one arrives at equations of motion for homogeneous fields 
0(t) and W(t) 



+ 3H^ + V' = j 



w+ I ij- 1 1 W 



W + 3HW 



f f 



2H' -H-~ 



f f 



W = 0. 



(4.1) 



(4.2) 



We assume / oc so that the last term of Eq. ( |4.2|) vanishes and there is a 
solution with W = 0. The other solution is decaying, W oc a~^, and corresponds to 
Bi = constant. With the action ( p.l| ) this can be set to zero by a gauge transformation 
and has no physical effect. We therefore take W(t) to be constant. (For hybrid 
inflation with the waterfall coupling to W both solutions are physical [01 but we 



*This result was first found in p6i by a different metliod. 
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assume that the decaying solution is anyhow neghgible.) It would make no essential 
difference if we allow / oc a" with a slightly different from —2, or more generally if 
we just assume 

i/iy ~- ( 2/f2-if^-^j W, if-^|W| < = |W|. (4.3) 



We assume that the right hand side of Eq. ( [4.1|) is small enough that the ap- 
proximation 3-ff0 ~ —V still applies. Then /'// ~ 2/\/2eM-p and treating the right 
hand side as a first order perturbation we have 

1 o i:t2tj/2 

3if0 ~ -V + —= -V. (4.4) 



The energy density of the gauge field is ||20| Pvk(x, t) = —f F^yF'^'' /A. Smoothed 
on a super-horizon scale this gives 

pw = ^H^W. (4.5) 

Using this result for the unperturbed energy density, we see that our assumption 
that the right hand side of Eq. ( f4.4| ) is dominated by the first term is consistent if 

^ - ^ « 1. (4.6) 

ep eMp 

which we assume. 



5. Perturbed universe 



Evaluating Eqs. (|2.5| ) and ( p. 61) to first order in 5(j) and dropping the right hand side 
of Eq. ( p^ we get 



+ 3H6<j) - V^6(f) = -V"t 



2eMp 



-5 



+ 3i7W^ 



( Wx + 3iJWx) 4- 



(5W + 3if5W - V^5W 



2W. 



5<p - 3H5(p - V^5<p 



V2eAfp 

Wx = W(x, t)=W + (5W(x, t), 



(5.1) 

(5.2) 
(5.3) 



where W is the unperturbed value. 

This ignores the metric perturbation (back-reaction). Since slow- roll inflation 
is supposed to be a good approximation, the effect of the metric perturbation in 
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Eqs. (|2.2|) and (2^) is significant only well after horizon exit, and is then given by 
Eq. (|3.8|) . Including it, Eqs. (|5.1|) and (p.2|) become well after horizon exit 



+ 3H6^ = - {V" + 6Hh) (50 + 
2 



+ 



(5Wx + 3H6W^ 



2eMp 
2W. 



+ 3HW^ 



2eMp 



2eMp 



- 3H6<f) + -He6(j) 



(5.5) 



where in Eq. ( p.4| ) the back- reaction —6H'^e6(j) is the same as in Eq. ( |3.17|) . 

Let us first set the right hand sides of Eqs. (|5.1| ) and (|5^) to zero. We saw in 
Section how the vacuum fluctuation of 6(f) is then converted at horizon exit to a 
nearly gaussian classical perturbation, and the same thing happens to the vacuum 
fluctuation of Wk Its left- and right-handed components have the same mode 
function W{k,t) = (j)o{k,t) in Eq. ( p.l5|) , which gives well after horizon exit 



{SW^SWl) = (27r)3 (6'^ - k'kn 6^{k + k')(27rVP)(if/27r) 



(5.6) 



Now we consider the effect of the right hand sides of Eqs. (|5.1|) and ( p.2| ). In 
the regime k ^ aH, spacetime curvature is negligible and we deal with field the- 
ory in fiat spacetime. The theory involves massless gauge bosons and the nearly 
massless infiaton particles {\V"\ ^ {k/aY' which propagate as nearly free particles 
(perturbative regime), justifying the initial condition (|3.9| ). 

To discuss the subsequent evolution, let us take (50(x, t) and 5W(x, t) to include 
only modes with k in some small interval, so that there is a well-defined epoch of 
horizon exit. In the (still quantum) regime k ~ aH, H is the only relevant scale and 
we have typical magnitudes 



(5.7) 
(5.8) 



Using these with Eq. ( ^.6| ), we see that the right hand sides of Eqs. ( |5.1| ) and ( |5.2| ) 
are much smaller than H^, and hence have only a small effect.^ It is therefore 
reasonable to assume that the perturbations become classical soon after horizon 
exit, with Eqs. (|3.16|) and ( |5.6|) initially a good approximation giving typical values 
|50| ~ \5W\ ~ H. 



^In contrast with the case for scalar field perturbations [|l5[, the non-gaussianity of 5W has yet 
to be evaluated, but we will assume that it still has a negligible effect. 



^The findings of this and the previous paragraph remain valid when the last term of Eq. (5.2) 
is included, and that term is negligible in the super- horizon regime that we are about to discuss. 
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The evolution of classical perturbations is described by Eqs. ( p.4|) and (|5.5| ) with 
k = 0. One of their solutions is 

3H5^ ~ - {V" + QHh) 54) + 5 (l^^(x, t)) (5.9) 



1 Q rr2 

= - {V" + ei/^e) + (2W ■ 5W + (5W ■ 5W) (5.10) 



2W 

5W ^ -^=— 5</>, (5.11) 



which implies \5(f)\ <^ -f^l^^l and |5W| ^ if|5W|. The self-consistency of this 
solution can be checked by inserting it into the right hand sides of Eqs. ( |5.1CI| ) and 
( ^.11|) . It is presumably picked out by the initial condition, just as in the case of 
slow-roll inflation. 

Since all of the quantities appearing in the second term of Eq. ( |5.1CI| ) vary slowly 
on the Hubble timescale, we will take them all to be constant, giving 



50k(W) = 50^(te„d) + -^1^ [m'))^ , (5.12) 



where 50o(iend) is the slow-roll result. 



6. Contribution of to the curvature perturbation 



Now we calculate C at the end of slow-roll inflation, using Eq. (|3.20|) . Smoothed on 
a super-horizon scale the energy density is p = p,^ + pw 



5p(x,tend) ^ 1^'50(X, tend) + \hH{W^) 



This gives C,]i{t 



end ) 



- V'6(po 

~ V^'(50o(x,tend) + im{k)HH{W^). 

Ck + Cif > where 

1 
2 

C{k) 
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x,tend) + lSN{k)H''5{W'') + ■^-H''5{W'' 



/■W 



'-Cik)5iW\ 



with 



W-5Wk 

QN{k) 



C{k) 



eM2 



(6.1) 

(6.2) 
(6.3) 



(6.4) 



(6.5) 



We can now calculate the spectrum and bispectrum of (w using Eq. (|5.6|) and 
the analogues of Eqs. ( |3.12| ) and ( p.l3|) . The result depends on the value of the 
unperturbed field W, and in particular on its magnitude W. 



-9- 



6.1 Working in a finite box 

An unperturbed quantity is the zero mode of its Fourier expansion, and at this point 
we need to remember that within the cosmological context that expansion has to be 
done within some box of finite coordinate size L Assuming nearly exponential 
inflation, the box leaves the horizon Nilk) = \n{kL) e-folds before the scale k. To 
avoid making unverifiable assumptions about an era that will never be constrained by 
observation (and, as we will see, also in to simplify the calculation) one should take 
koL to be as small as is consistent with the requirement that the periodic boundary 
condition implied by the use of the Fourier expansion have a negligible effect (minimal 
box). Demanding say 1% accuracy in the calculation, it should be enough to choose 
koL ~ 100 corresponding to ln{koL) ~ 5. After choosing L one writes for a given 
quantity 5'(x, t) = g(t) + 5g{'x.,t), where the 'unperturbed' value g{t) is the average 
of g within the box. 

Expectation values like ( p.lO|) are in general defined with respect to an ensemble 



of universes, one of which is the observed universe. But under the usual assumption 
that perturbations originate as a vacuum fluctuation the translation invariance of 
the vacuum makes them translation invariant. As a result the expectation values 
can be defined as spatial averages within a single realisation of the ensemble. Thus 
the {{SgY) can be defined as the spatial average within the box, {Sgi^ + X.)6g{'x)) 



can be defined as the average with respect to X and so on ||13 
Keeping only classical modes, Eqs. ( p.l6| ) and (^]6|) give 



k 

2 



{{S<P))) = J V4k)dk/k = \n{koL){H/27Ty (6.6) 

(|(5W|2) ^ 2\n{koL){H/27T)\ (6.7) 

For the minimal box size this corresponds to typical values ~ |(5W| ~ H. 

6.2 The value of W 



As was discussed in ||21|, there are two possible viewpoints about the magnitude of an 
unperturbed field like W.'' One is to regard the magnitude as a free parameter, 
that one gets to choose just like the masses and couplings appearing in the action. 
The other viewpoint is to estimate the likely value of W, assuming that we are at a 
typical location within a box whose size M is very much bigger than the size L of 
the minimal box within which the calculations are done. 

Adopting the second viewpoint one has to make an assumption about the evo- 
lution of the universe long before the observable universe leaves the horizon. The 



^These apply to any non-inflaton field that acquires a perturbation from its vacuum fluctuation. 
The inflaton field is an exception because the inflation model and the cosmology determines its 
value N{ko) before the end of inflation. We ignore the issue of anthropic selection, assuming that 



statistical anisotropy (like non-gaussianity 13 ) is neither favoured or disfavoured in that respect. 
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usual assumption is that there is almost exponential inflation, beginning Nm e-folds 
before the observable universe leaves the horizon with Nm fairly large. One can then 
estimate the likely value with suitable assumptions about the relevant physics during 
those e-folds. In our case, let us assume that the dependence / oc continues to 
(at least approximately) hold during those e-folds. Then, after smoothing W on the 
scale L Eq. (|5.6|) gives for the average within the exponentially inflated patch 



(|5Wn ~ 2\Yi{M/L){H/2Txf ~ 2Nl{k){H/2nf , (6.8) 

The expected value of if our location is typical and therefore 

= WIj + \n{M / L){H /2-Kf > \n{M / L){H /2-Kf = A^A/(A;o)(^/27^)^ (6.9) 

where Wm is the average within the inflated patch. We conclude that if we occupy 
a typical location within an inflated patch that left the horizon many e-folds before 
the observable universe, then W ^ H. 

6.3 The case W :^ H 

If W ^ H, the second term in the square bracket of Eq. (|6.4| ) can be treated as a 
first-order perturbation. Then, assuming that (w gives the only contribution to the 
anisotropy of V(, Eq. (39) of P| gives 



Pc(k)=n(^) 

with 



l + g,{k) (w-k 



(6.10) 



9*ik) = -VUk)/V^ik) (6.11) 
- ^^(^)^7^V (6.12) 



Here, is given by Eq. ( ^.26] ), and is independent of k since we are taking H and 
e to have negligible time-dependence. 

The anisotropy is of the form (^TTj) but with a strongly scale-dependent g*{k) oc 
N{k). This has not been compared with observation but the constraint is presum- 
ably similar to the scale-independent case, currently \g^:\ < 10~^. Another con- 
straint comes from the strong scale-dependence of Vc^^ik), giving a contribution 
n{k) D —g*{k) to the spectral index. Using Eq. ( p.23|) , this requires barring a can- 



cellation |5'*(/co)| ^ 0.03. In any case, one certainly needs |5'*(fco)| ^ 1, which with 
Eq. ( p.l3p is stronger than Eq. ( [4. 61 ), justifying the latter. Our assumption that 



the anisotropy of the expansion has a negligible effect is also justified, because that 
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anisotropy is presumably only of order {pw/p), which presumably gives a contribu- 
tion {pw / p){'Pc4,/'Pci^)) 9* which is smaller than the one that we have calculated. 
Using V^iko) = (5 x 10-^)2, 



2 /Tjr\ 2 



If we were to assume "P^^ = V( as in [|], the observational bound < 10^^ leads 
to two conclusions as those authors notice. First, our assumption W ^ H would be 
only marginally allowed. Also, from Eq. ( p.9|) , one sees that Nm to be very large if 
we occupy a typical location. 

Now we calculate the contribution of (w to the bispectrum on cosmological 
scales. To simplify the calculation we set N{k) = N{kQ). Taking VQ{k) to be scale- 
independent and defining 

^/NL(kl,k„k3)^^^M, (6.15) 

Eq. (41) of H gives 

/NL(ki, k2, kg) = /nl (1 + /ani(kl, ks, kg)) , (6.16) 

where the constant prefactor is given by 

f 5 gl^ko) 

= -im{kMko){VcjVc), (6.18) 

and 

= -(A-kO^-(A.fa)^+(k..k.)(A.k.)(A.M ^ ^ ^^^^^^ 

This expression has not yet been confronted with data, but the bound is presumably 
similar to the one that takes g^:{k) to be constant, currently |/nl| ^ 10^. 

6.4 The case W <H 

Consider now the case IV = 0. Since there is no preferred direction, is statistical 
isotropic. Using Eq. (40) of p| we have 



' Qw _ Qw 
Vc ~ Vc 



\c\k)(^^\n{kL)/V^ (6.20) 



C ' c 

2 



?>N^{k)V^\Yi{kL)(^^^^^ <1. (6.21) 
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Using Eq. (5^) and the analogues of Eqs. ( p.l2 ) and (|3.13 ) we find 



Jnl — /, 



loop 
NL 



2 fHyinjkL) 

3 V27r; {V^y 



1 + 



V 



kf/k 



3/U3 



3^/2 fn 



4 V V. 



3/2 



C 



9iV3p^ ) HkL) 



1 + 



1 + 



ki-k, 



V 



< 1 



(6.22) 



(6.23) 



(6.24) 



In the general case, Vc^^ and /nl are the sum of the 'tree' and 'loop' contributions. 
6.5 How to make the calculation more accurate 

The calculation that we have presented is sufficiently accurate, unless and until 
statistical anisotropy is observed. If that does happen a more accurate calculation 
may be required. 

Such a calculation should in principle begin with the generation of the classi- 
cal perturbations 50 and (5W from the vacuum, along the lines of It may in 
principle generate significant time-dependence and correlation for these quantities, 
when they first become classical. But the analysis of Section ^ suggests that it will 
instead confirm the result obtained there; nearly time-independent and uncorrelated 
perturbations with the spectra (if/27r)^ defined in Eqs. (p.lO|) , ( p.l6|) , and ( |5.6|) . 

After the perturbations become classical, their evolution is given by Eqs. (|5.4] ) 
and ( p.Sp , which can be solved numerical to determine the correlators of the pertur- 
bations at the end of infiation. That allows one to determine 5p, and working to first 
order in C, we can use Eq. ( |3.2CI| ) to determine and /nl- 

Instead of working with the perturbations, one can use the 5N formula |^ 

C(x,t) = C0(x,t) + Ciy(x,t)+ + \N^'t>ii)\^<P*i?^)f + ^iV^.50*(x)5W(^^25) 

i 

C^i^,t) = N^it)5M^), (6.26) 
Cwi^,t) ^ ^iV,(t)5H^*(x) + i5^iV,,(t)51^*(x)5iy,*(x). (6.27) 



Here N is defined by Eq. (|3.19| ) with the epoch of horizon exit for the scale of 
interest, and the subscripts on denote partial derivatives evaluated at the unper- 
turbed point in the field space. But as we now argue, the perturbative approach of 
the previous paragraph is expected to be adequate, i.e. the first-order Eq. (|3.2CI|) is 
expected to be adequate. 
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The validity of Eq. ( p.20|) for C,^ is a standard result but we need to justify its use 
for C,w The second order correction is presumably of order Qy, corresponding to a 
tiny fractional correction of order C,y/ that can certainly be ignored for the evaluation 
of Vc^yy. To see whether it can be ignored for /nl, let us first pretend that W is a 
scalar field. Then, setting C{k) to a constant and assuming that the first term of 
Eq. ( |6.4| ) dominates, Eqs. ( |6.4] ) and ( |6.17D give^ 



C(x) = Cg(x) + ^/nl ( ) C(x), (6.28) 



3 . / n W2. 



where Cg = CW5W. The first-order formula will therefore be adequate unless 

/NL<(^y. (6.29) 

Since we need |/nl| ^ 1 for it to be observable, we conclude that the first-order 
formula will be adequate for the evaluation of /nl unless |/nl| ~ 1 o.nd (w is the 
dominant contribution to (. Except for the second proviso this is a standard result, 
that was first recognised in the context of the curvaton scenario |2^, . 
Keeping the vector nature of W, Eqs. ( |6.4| ) and ( |6.17D give 



C.(x) = Cg(x) + ^/nl C|(x) (6.30) 

Cg(x) =CW-5W(x)) (6.31) 
Cg(x) = Ciy|5W(k)| (6.32) 

At a typical location, ~ Q, and the previous conclusion about the validity of the 
first-order formula still applies. 



7. Conclusion 



Working exclusive with the classical perturbation of W, we have presented a rather 
complete calculation of the contribution to (w that is generated during slow-roll 
inflation. In the regime W ^ H we have reproduced the result of |]10|, ^ for the 



spectrum (Eq. ( |6.13D ) and of P] for the bispectrum (Eqs. ( |6.16D and (|6.18| ), discussing 
for the first time the assumptions that are needed to obtain it. 

The field W(x, t) may survive after the end of slow-roll inflation, and generate 
further contributions to (. The effect of these is to change C{k, t), so that Eqs. ( |6.11| ) 
and (|6T?|) as well as Eqs. ( lOOD , ( jO^ ), and ( IQgD still hold but with C the final 
value. 



*If instead the second term dominates, it is solely responsible for /nl and the first-order formula 
will certainly be adequate. 
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A contribution (end may be generated during the waterfall that ends hybrid 
inflation. It can be calculated from the 'end-of-infiation' formula, provided that the 
waterfall is sufficiently brief which requires H < 10~^Mp('P^^^^/'P^)^/^.^ The result 
is 

'^cnd = 1=— (7.1) 

y/IeMping 

If the gauge symmetry is spontaneously broken after inflation, the curvaton 
mechanism may generate another contribution after inflation. With the simplest 
assumption that W is nearly time-independent until it begins to oscillate 

where 'dec' denotes the epoch just before the decay of W. 
Including both these contributions we have 

C = Csr + C'end + C'curv; (7-3) 

where Csr = QN{k)/eMp is the contribution generated during slow-roll inflation. The 
extra contributions allow the non-gaussianity to be observable even if W H. 



We close with an important comment, stemming from Eq. (|6.14|) . On the rea- 
sonable assumption that there were some large number Nm of e-folds of inflation 
before the observable Universe left the horizon, one expects W/H ~ ^ 1- But 
then Eq. (|6.14|) is compatible with the observational bound on ^f* only if V^^ ^ "P^; 
in other words, if the the observed curvature perturbation is mostly generated after 
slow-roll inflation by for instance the end-of-inflation or curvaton mechanism. Our re- 
sults therefore suggest a dichotomy regarding the generation of observable statistical 
anisotropy of the curvature perturbation. Either the entire curvature perturbation 
( (both the isotropic and the statistically anisotropic part) is likely to be generated 
during slow-roll inflation, or it is likely to be generated afterwards. 
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